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Abstract 

Initial data for the spherically symmetric Einstein- Vlasov system is 
constructed whose past evolution is regular and whose future evolution 
contains a black hole. This is the first example of initial data with these 
properties for the Einstein-matter system with a "realistic" matter 
model. One consequence of the result is that there exists a class of 
initial data for which the ratio of the Hawking mass m = nir) and 
the area radius r is arbitrarily small everywhere, such that a black 
hole forms in the evolution. This result is in a sense analogous to the 
result (T5] for a scalar field. Another consequence is that there exist 
black hole initial data such that the solutions exist for all Schwarzschild 
time t £ (—00, 00). 



1 Introduction 



An important question in the study of gravitational collapse is to identify 
physically admissible initial data, and it is natural to require that the past 
evolution of the data is regular. However, most of the existing mathematical 
results which ensure a regular past also ensure a regular future, cf . \22\ [T6] , 
which rules out the study of the formation of black holes. The exceptions 
being the classical result for dust [19] , where some classes of solutions have a 
regular past, and the recent result [T7] for a scalar field. In the latter work, 
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which in part rests on the studies |14[ I15|. initial data whose past evolution is 
regular and whose future evolution forms a black hole is constructed. Now, 
neither dust nor a scalar field are realistic matter models in the sense that 
they are used by astrophysicists. Dust is a perfect fluid where the pressure is 
assumed to be zero, and a scalar field is merely a toy model. Thus, there is 
so far no example of a solution to the Einstein-matter system for a realistic 
matter model possessing a regular past and a singular future. 

In this work we consider collisionless matter governed by the Vlasov 
equation, cf. P] and [25] for an introduction. Although this is a simple 
matter model, it has rich dynamics and many features that are desirable of 
a realistic matter model. Indeed, it allows for anisotropic pressure, there 
is a large number of stable and unstable spherically symmetric and axi- 
ally symmetric stationary solutions [8j [21], there is numerical support that 
time periodic solutions exist [10], it behaves as Type I matter in critical 
collapse pHl HS1 [24] , and it is used by astrophysicists [12] . The following 
theorem is the main result of the present paper. 

Theorem 1 There exists a class of initial data J for the spherically sym- 
metric Einstein- Vlasov system with the property that black holes form in 
the future time direction and in the past time direction spacetime is causally 
geodesically complete. 

We will see that a consequence of this result is that for any e > 0, initial data 
can be constructed with the property that the ratio m/r of the initial Hawk- 
ing mass m= m(r), and the area radius r, is less than e everywhere, such 
that a black hole forms in the evolution. We formulate this as a corollary. 

Corollary 1 Given e > 0, there exists a class J r of initial data for the 
spherically symmetric Einstein-Vlasov system which satisfy 



for which black holes form in the evolution. 

This result improves the main result of [6] and is analogous to the result [15] 
in the case of a scalar field where conditions on the data are given which 
ensure the formation of black holes. These conditions give no lower bound 
on 2m/ r but involve other restrictions. Another consequence of our result is 
the following corollary. 

Corollary 2 There exists a class J s of black hole initial data for the spher- 
ically symmetric Einstein-Vlasov system such that the corresponding solu- 
tions exist for all Schwarzschild time t £ (—00,00). 
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In the future time direction this corollary was shown in [6] , the improvement 
here is that the solutions exist on the entire real line. 

The present result relies in part on the previous studies [2], [5] and [6], 
which now will be reviewed. In [5] global existence in a maximal time gauge 
is shown for a particular class of initial data where the particles are moving 
rapidly outwards. One of the restrictions imposed on the initial data is that 



where the constant ko is roughly 1/10. The situation considered in [B] is in a 
sense the reverse since the initial data is such that the particles move rapidly 
inwards and the quantity sup r 2m/ r is required to be close to one. The 
main result of [6] is that data of this kind guarantee the formation of black 
holes in the evolution. The analysis in [B] is carried out in Schwarzschild 
coordinates, i.e. in a polar time gauge. Now, particles that move inward in 
the future time direction move outward in the past time direction. It is thus 
natural to try to combine these two results with the goal of constructing 
solutions with a regular past and a singular future. The conditions on the 
ratio 2m/ r are clearly very different in [5] compared to [6], and moreover, 
the Cauchy hypersurfaces are different since a maximal time gauge and a 
polar time gauge are imposed in the respective cases. The main reason why 
a maximal time gauge is used in [5] is due to the difficulties related to the 
so called pointwise terms in the characteristic equations in Schwarzschild 
coordinates. In [2] the problem of global existence for general initial data is 
investigated under conditional assumptions on the solutions. The analysis 
along characteristics is applied to a modified quantity for which the problems 
with the pointwise terms in Schwarzschild coordinates do not appear. 

In the present work we combine the strategies in [2] and [5] and show 
global existence for rapidly outgoing particles in Schwarzschild coordinates. 
In particular the result in [5] is improved by showing that the restriction (jl.ip 
can be relaxed, and for sufficiently fast moving particles 2m/ r is allowed to 
be arbitrarily close to one. By combining this result with the result in [6| 
we are then able to construct data whose past is regular and whose future 
contains a black hole. 

The outline of the paper is as follows. In the next section the spheri- 
cally symmetric Einstein- Vlasov system is introduced. Global existence for 
rapidly outgoing particles is shown in section 3 for two different sets of initial 
data which are adapted to Corollary 1 and Corollary 2 respectively. Finally, 
in section 4 the proofs of Theorem [JJ Corollary 1 and Corollary 2 are given. 
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2 The Einstein-Vlasov system 

For an introduction to the Einstein-Vlasov system and kinetic theory we 
refer to pQ and [25] , and for a careful derivation of the system given below 
we refer to [20]. In Schwarzschild coordinates the spherically symmetric 
metric takes the form 

ds 2 = _ e 2Kt,r) df 2 + e 2X(t,r) dr 2 + ^^2 + ^2 ^ 

The Einstein equations read 

er 2X {2r\ r - 1) + 1 = 8irr 2 p, (2.2) 
e- 2X (2rfi r + 1) - 1 = 8irr 2 p, (2.3) 
\ t = -4vrre A+ ^j, (2.4) 

e- 2X (fi rr + (fi r - X r )(fi r + -)) - e- 2 ^(A tt + X t (Xt - fit)) = 8tt Pt . (2.5) 

The indices t and r denote partial derivatives. The Vlasov equation for the 
density function / = f(t, r, w, L) is given by 

8 t f + e^ x ^d r f - (Xtw + er x HrE - e^ x ^)d w f = 0, (2.6) 

where 



E = E(r, w, L) = y/1 + w 2 + L/r 2 . (2.7) 

Here w E (— oo, oo) can be thought of as the radial component of the mo- 
mentum variables, and L £ [0, oo) is the square of the angular momentum. 
The matter quantities are defined by 

p (t,r) = —l I Ef(t,r,w,L) dwdL, (2.8) 

r J-oo JO 

[oo /'OO w 2 

P^ r ) = ^i / - W f(t,r,w,L)dwdL, (2.9) 

J — oo J 

,-00 /-00 

j(t,r) = - / wf(t,r,w,L) dwdL, (2.10) 

r J —oo J 

/"oo /"oo r 

PT(t,r) = — [ / -f(t,r,w,L)dwdL. (2.11) 
J-oo Jo tj 

Here p,p,j and are the energy density, the radial pressure, the current 
and the tangential pressure respectively. The following boundary conditions 
are imposed to ensure asymptotic flatness 

lim X(t,r) = lim n(t,r) = 0, (2.12) 
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and a regular centre requires that 



A(t,0) = 0. (2.13) 

We point out that the Einstein equations are not independent and that e.g. 
the equations ([13]) and (J23J) follow by ([22]), ([23]) and ([2"U]) . 

o 

As initial data it is sufficient to prescribe a density function / = 
f(r, w, L) > such that 

J ^r ] 2 p{j 1 )dr 1 < r -. (2.14) 
Here we denote by p the energy density induced by the initial distribution 

o 

function /. This condition ensures that no trapped surfaces are present 

o 

initially. Given /, equations (|2.2p and (|2.3h can be solved to give A and 

o 

p, at t = 0. We will only consider initial data such that f(r,-,L) = if 
r < e, for some e > 0, or if L > L + , for some L + > 0. If the initial data is 
C 1 ([e, oof, ] — oo, oof, [0, oof) we say that it is regular. 

Let us write down a couple of facts about the system (|2.2p - (|2.12p . A 
solution to the Vlasov equation can be written 

f(t, r, w, L) = f(R(0, t, r, w, L), W(0, t, r, w, L),L), (2.15) 

where R and W are solutions of the characteristic system 

^ _ >-A)(s,fl) W (2 16) 

ds~ E(R,W,LY { ' 

^ = -X t (s,R)W -e^- x)is ' R) p r (s,R)E(R,W,L) 
as 

■ (u-X)(sM) L (9 1 7^1 

+C R S E(R, W, L) ' {ZAn 

such that (R(s,t,r,w,L),W(s,t,r,w,L),L) = (r,w,L) when s = t. This 

o 

representation shows that / is nonnegative for all t > 0, ||/||oo = ll/lloo, and 
that f(t, r, w, L) = if L > L + . The Hawking mass m of the sphere of area 
radius r is given by 

m(t,r) = 4m- rf p(t,r])dr], (2.18) 
J o 

and by integrating (|2.2p we find 

= 1 - 2m(t ' r) . (2.19) 
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A fact that we will need is that 

fi + A < 0. 

This is easily obtained by adding the equations (|2.2[) and (|2.3p . which gives 

A r + /i r > 0, 

and then using the boundary conditions f)2. 12|) . Furthermore, from ()2. 19[) 
we get that A > 0, and it follows that fi < 0. We also introduce the notations 
p, and p. From equation (|2.3j) and f|2. 12|) we have 

f°° m(t n) f°° 
(j,(t,r) = - v ' " e 2X drj - / 47rr/pe 2A drj =: fx + p.. (2.20) 

Jr V Jr 

We will need an expression for fit. By (j2.4j) and (j2. 19|) it follows that 
mt(t,r) = —A7rr 2 j(t,r)e fJ- ~ x , and we obtain 



p t {t,r)= / 47ri(t,r/)e^ +A )^e 2A ^dr ? . (2.21) 

Jr 

An important quantity is the ADM mass M, given by 

/•OD 

M = 4ir r 2 p(t,r)dr. (2.22) 
Jo 

The fact that it is conserved follows by using (j2.4|) and (j2. 19|) . 

The following result is given in [5] but since the proof is short we include 
it for completeness. By a regular solution we mean a solution which is 
launched by regular initial data with compact support. 

Lemma 1 Let (/, [/,, A) be a regular solution to the Einstein- Vlasov system. 
Then 

POO 

/ 4Trr(p + p)e 2X e^ +x dr < 1, (2.23) 
Jo 
f°° m 

/ (-J + 4Trrp)e 2X e fl dr < 1. (2.24) 
Jo r 

Proof. Using the boundary condition (|2.12p we get 

poo J 



1 > l-e^- A (t,0) = / —e^ +x dr 
' dr 

(fir + X r )e^ +X dr. 

The right hand side equals (|2.23[) by equations (|2.2p and (|2.3p which com- 
pletes the first part of the lemma. The second part follows by considering 
e M instead of e M+A . 
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Finally, we note that in [22] and [2] local existence theorems are proved for 
compactly supported and non-compactly supported initial data respectively, 
and it will be used below that solutions exist on some time interval [0, T[. 



3 Global existence for outgoing matter 

In order to understand the mechanism behind the global existence result for 
outgoing matter we recall the example in [5] and consider the much simpler 
Vlasov-Poisson system which is the Newtonian limit of the Einstein- Vlasov 
system. Due to the spherical symmetry the maximal force experienced by a 
particle at distance r from the origin is — M/r 2 in the Vlasov-Poisson case, 
where M > is the total mass of the particle ensemble. Hence along any 
particle trajectory 

d (\ 2 M\ M ( M 

~r ^w = ww H — = w w H — T > 



dt \2 r J r 2 \ r 2 

as long as its radial velocity r = w = x- v/r>0. Hence 

1 . M 1 9 . , M 
-W 2 (t) --r-r> -W 2 (0) 



2 v ' r(t) ~ 2 w r(0) 
and 

-w 2 (t) > -w 2 (0) 



2 w - 2 w r(0) 

on any time interval on which w(t) remains non-negative. Now let w\ > 
be an initial lower bound for the radial velocities of the particles in the 
ensemble, n > an initial lower bound for their distance from the origin, 
and assume that 

1 n M 

Wi := -w{ > 0. 

2 n 

Then as long as a particle is moving outward, 

w(t) > W u r(t) >n + Wit. 

This implies that all the particles keep moving outward for all future time. 

Let us turn back to the spherically symmetric Einstein- Vlasov system. 
The notation below is adapted to the notation in [6] since the aim is to show 
that the initial data we construct overlap with the initial data in [6J. Two 
different sets of initial data, adapted to Corollary 1 and 2, will be considered 
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and two similar results on global existence will be shown below; Theorem 2 
and Theorem 3. Let us point out that if the only goal had been to improve 
the global existence result for rapidly outgoing particles in [5] then we could 
have considered a simpler class of initial data analogous to the data in [5], 
cf. Remark 2 below. 

o 

Let < ro < r\ be given and put M = r\/2. Let f s be data of a steady 
state supported in [rg,ro] and let 

pro 

M in := / iTrr 2 p s (r)dr, (3.25) 

H 

o 

where p s is induced by f s . The results in [1] guarantee that such steady 
states exist if r' is sufficiently small, and moreover that 

2mir) 8 

sup — — < -. (3.26) 

0<r<r T 9 

This implies in particular that 2M; n /ro < 8/9 so that M > 9M; n /8. Let 
M out := M — M in . Let R x > n be such that 

fli-ri<^=^, (3.27) 
b 

and define 

i?o := -(n + Ri). 
Let L + > and let W* > be such that 

|W*|>l + ^±. (3.28) 
Ho 

Let W- > satisfy 



\W-\ e 2Ro(1 ~^ } (1 - ^£)3/2 > 3 |^|. (3.29) 
ito 

We remark that since W* > and > in this section, the modulus is 
superfluous but it will be needed below. The same remark applies to the 
time variable which is non-negative in this section but which will be non- 
positive below and we therefore in some situations write the modulus of the 
time variable, cf. ()3.32p . We are now in a position to specify the initial 

o o o 

data. Let f=f s + fm be initial data of ADM mass M such that 



supp/ m C [R ,Ri] x [W-,oo[x[0,L 



+J > 
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and 

rRi 

4nr 2 $ m {r)dr = M ont , (3.30) 

iRo 



f 



where pm is the induced energy density by f m . 

Remark 1: Note that the condition (j4.4p can be arranged by first choos- 
ing h m such that supp/i m C [Ro,R\] x [W-, oo[x [0, L + ] and then choosing 
an amplitude A £ M + such that f m := Ah m satisfies (|4.4p . 

Before stating the main result in this section we define 

IWJ " 2M y^W\ (3.31) 



y/1 + + L+/R* Ro 



o 

Theorem 2 Assume that r , ro, fi, M{ D , M, L + , Rq, R\, W*, W- and f are 
given as above, and consider a solution f of the system i2.2\) - [27h]) . launched 

o 

by f, on its maximal existence interval [0,T[. Then T = oo, and 

supp/ m (£) C [.Ro + |t«*|,oo[x[W*,oo[x[0,L+], (3.32) 
and the resulting spacetime is future causally geodesically complete. 

Remark 2: The initial data in the theorem is adapted to match the initial 
data constructed in [6]. However, it is important to note that the presence 
of the steady state given by f s is not needed if the only aim is to construct 
initial data for proving global existence and geodesic completeness. Hence, 
by taking f s = and disregard the parameters r' ,r$ and M m and thus 
consider a simpler class of initial, then Theorem 2 can be directly compared 
with the result in [5], and it can be seen to be an improvement of this result. 

Proof: We first notice by following the arguments in [6], that the only 
way the matter in the outer region r > Rq can affect the static solution is 
via the metric function \i. By dropping the time derivatives in the Vlasov 
equation we see that in the remaining equation the factor e^ _A can be can- 
celed. Hence, the static Einstein- Vlasov system is formulated in terms of 
the quantities /, A and [i r and not [i itself. Therefore /, A and [i r remain 
time independent for r < r§. The arguments of the proof will therefore 
mainly concern the outer matter given by f m . In particular we will see 
that the outer matter which initially is moving outwards will continue to 
move outwards, and therefore there is no direct interaction with the steady 
state. However, in the last part of the proof which concerns causal geodesic 
completeness the steady state will have an influence. 
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Let [0, t\ [ be the maximal time interval such that for t G [0, t\ [ and 
(r,w,L) G supp/ m (t), w > W*. By continuity t\ > 0. Suppose that 
t\ g]0,T[, then we must have u; = W* for some w G supp/ m (ii), but we 
will show that w > W* for all w G supp/ m (ii). Thus t\ = T and since the 
matter stays strictly away from r = it follows that T = oo in view of [2] 
or [23]. 

Consider a characteristic (R(s),W(s), L) with i2(0) G [Rq,R\] and let 
G(i) := E(R(t),W(t),L) + W(t). 
Note that G > 0. We have by (^TUD and ([2TT7D that 



dG(s) 
ds 



X t (s,R(s))- 



W(s) 



E(R(s),W(s),L) 

Le (V-\)(s,R(s)) 



+ /i r (s, J R(s))e^- A » s '^ s » G{s) 

+ R{s) 3 E{R(s),W(s),L)' (3 ' 33) 

Below we will often suppress the arguments but it should be clear that 
R = R(s), n r = n r (s,R(s)) etc. The quantity G, which was first introduced 
in [2], is not suitable for the purpose here and the quantity that we will 
consider is G{t)e^' R ^ (1 - 2M/R(t)). We have 



d_ 

ds 



(Ge^l 



2M 
~R~ 



)) 



' W 



« 3 £ 
2MW 

~R?~E 



2M 



Ge^l 



) + [Ar 



— i 



2M 



2M S 
~R' 



e^- x Ge fl 



At— + fi r e^ x 



+^-e"(l 



2M, 



+ 



« 3 £ 
~R?~E 



e^ A Ge A . 



ft. G<"[1 -^-) 
(1 + Jr) m +A 



e^(l 



2M 



i? ' 
(3.34) 



Here we used that 
fi r —e»- x G 



■\i T e 



:i + 



■/?? 



E r "~ E R 2 

Consider the second last term in (|3.34|) . In view of (|2.19p we have 

2M 



1 



R(t) 



< e -2X(t,R(t)) 
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and we get 



IT 1 D2 V d ' — D2 



Now, since > W* on [0, ti], and since > Ro, we get in view of f)3.28|) 
that 

(I + W) ™ xa 2MW x a 

E IV R 2 E 

L 



, (1 + -m) 2WG\ m ,. \ n / \ 

> ( " + — ) ^e"" V > 0. (3.35) 



Hence we have 



This implies that 

2M 



G^(l-^p). (3.36) 



G(t 1 ) e A(*i--R(ti))(i 



> e 



-/ *i[A t ( s> fl( s ))f +Ar(s,R(s))e('*- A )C^ Ji ( s ))-At(s,ii( S ))]<i SG / ) e /i(O,i?(O))/ 1 



2M 
(3.37) 



Let 7 be the curve 

7 :={(*. r):0<*<ti,r = 
The time integral in (|3.37p can be written as 

j e (-n+\)(t,r) Xt ^ r ^ dr+ ^- A )(*' r )/i r (i,r)-/lt(t,r)) di. (3.38) 

We will apply Green's formula in the plane to this curve integral. Let 
Roo > Ri + H, so that /(-,i?oo,-) = for t 6 [0, ti]. Let T be the closed 
curve 

r = 7 + c tl + Coo + c , 

where 

C tl := {(i,r) :t = ti, R(h) <r< R^}, 
C :={(t,r):t = 0, R(0) < r < R^}, 



11 



and 

We have 



Coo ■= {(t,r) : h > t > 0, r = R QO }. 

j> e~^ +x X t dr + (e^fa - fa) dt 
8 t (e^ +x X t ) - d r (e^ x fi r - fa) drdt 



8 t (e-^Xt) - d r ( e"" V - At) drdt + d r ( e^ x 



m 



drdt. 
(3.39) 



By using (|2.5p and (|2.2ip this identity can be written 
j> e~^ +x X t dr + {e^~ x fi r - fa) dt 



2m 



4ir(p — p) — 8irpx — 47re 2A j ) drdt 



+ 



/,_N ' (fir + A ; ;— 7 • \7T/> — ) drdi 



n 

4vre^ +A 



(p + p)e 2A — +p-2p T -je 2A 
r 



drdt. 



(3.40) 



Here we used that \x T + A r = 47rr(/? + p)e 2X . Since u; > W* > on [0, ii] we 
have in view of (|2.9p and (|2.10p that j > p, and we have 



p(l + e 



2A 



— = pe 2A 1 + — < pe 2A . 



Hence, by dropping the term involving pt due to sign, we thus obtain 
j> e->* +x X t dr + (e>*- x fir- fa) dt< j J 4vr A pe 2X ™ drdt 



< 



Aire^pe 



m 



-drdt. 



(3.41) 



Here A = {(t, r) : < t < ti, i?o ^ r < an d the last inequality follows 
since the integrand is nonnegative and Q C A. Next we estimate p in terms 
of j. Since for w 6 supp/ m (i), < t < ti, we have u> > W* > 1 + y/L + /Ro, 
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we get for r > Rq, 

roo poo poo poo 

p(t,r) < — / / fdLdw-\ — g / / wfdLdw 
r J-oo Jo r J-oo Jo 

+ - 2 T f°°—fdLd W 



>' 1-oaJO 

"CO roo 



^_ poo poo 

<3— / / wfdLdw = 3 j(t,r). (3.42) 

r J-oo Jo 

We estimate the right hand side in ()3.41 j) by using the Vlasov equation from 
which it follows that 

| (r 2 e 2X p(t,r)) = (rV+*j) - re^ x 2je 2X ^. (3.43) 

Since j(t, Rq) = and j(t, Roo) = 0, this gives 

Roo fRoc 

r 2 e 2A(ti,r)^ ti>r ) dT ._ / r 2 e 2X ^p(0,r)dr 
Ro JRo 

re> M+x 2je 2X —drdt. 



A 



Thus, we get 

/ / re» +x 2je 2X -drdt < H r 2 e 2X p(0,r)dr < (3.44) 
J JA r Jr 47F 1 ~ Tfe 

where the last inequality follows in view of (|2.19p and (|2.22p . Using inequal- 
ity (|3.42p we therefore have 

eV +x pe 2X ^drdt<4ir J J e » + Hje 2X ™drdt 

< 4vr-^ ! f re» +x 3je 2X -drdt 
~ Ro J J A r 

< M 9Af . (3.45) 

-2^(1-^) 

In order to obtain an estimate for (|3.38p it remains to estimate the boundary 
terms since 



'7 Jr JCt 1 JCoo J Co 

First we notice that the curve integral along Coo vanishes since both p and 
j, which determine p r and fit, vanishes for r = R^. Since j > we have 
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that At < 0, which implies that the integral along Co can be dropped due 
to sign noticing the orientation of Cq. The term along Ct x can be estimated 
as follows. 



Ci 



< 



Roo 

R(ti) 
M 



,(-^+mi,r)x t (t ur )dr < / 4TTre 2X \j{h,r)\dr 



Ro 



r> / 1 2M n 



(3.46) 



In the first inequality we used (|2.4p and in the second we used (|2.19p . the 
fact that \j\ < p, and (12.22p . We have thus obtained the estimate 



^ +A A< dr + (e"- A Ar - fL t ) ds < 



5M 



7 



2*0(1-^)" 



Inserting this into the main inequality we get 

-5M 



G{h)e 



i 



2M 



> e 



»>G(0)e 



£(Q,-R(0)), 



1 



2M . 

W)' 



Noticing that fi is monotone in r and nonpositive, and that R(0) > Ro, we 
obtain the inequality 



—oivi 

G(ti) > e 2 " 0(1 -^ ) G(0)e^)(l - — ) 

Ro 



-5M 



Ro Ro 



Here we made use of the estimate 

"00 



2M 
i?0 ' 



V 



(3.47) 



(3.48) 



We have that G(0) > W(0) > W_, and in view of (I3T28I) we also have 
3W(t) > G(t) on [0,£i]. We now use the condition (pT29D and obtain 

3W(h) > G(h) > 3W*. 

Thus W[ti) > W*, and necessarily we have ii = T. As was pointed out in 
the beginning of the proof, since matter stay strictly away from the centre of 
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symmetry, T = oo, cf. [2] or |23| . Let us next show that ()3.32j) holds. From 
the characteristic equation (|2,16p we can conclude that R(t) > i?o + \ ft*t\ in 
view of the estimates 

P -x(tMt)) -fi_ 2m(t,/Z(t)) y/a / _ 2MU/2 

V 1 i?(t) J -V rJ ' 

and fj,(t,R(t)) = £t(t,R(t)) + p,{t,R(t)) > p,(t,R ) + p,(t,R ) where fi(t,R ) 
is estimated by (|3.48p and 

£fo Ho) = - / 4ir VP e 2X ^ d V > - 

jRn -Ko(J 



The latter inequality is analogous to (|3.46p . We remark that the estimate 
R(t) > i?o + is rough and can be improved by using arguments from [9]. 
In order to complete the proof of the theorem we have to show that any 
causal geodesic is complete. We follow the argument in [5] and introduce 
the coordinates 

x°=t, x 1 = r sin 6 cos (f>, x 2 = r sin 8 sin <fi, x 3 = rcos9. 

In these coordinates the metric becomes 

5oo = -e 2A \ goa = 0, g a b = &ab + (e 2A - 

where Latin indices a, b run from 1 to 3 and x a = S ao x b . Let us now con- 
sider an arbitrary future directed, time- like or null geodesic, i.e., a solution 
(x 7 (t),p 7 (t)) of the geodesic equations 

dx' 7 
dp 1 



where Greek indices 7, 5, e run from to 3, Tj e are the Christoffel symbols, 
and 

p° > 0, g 7 ipV = < 0. 

We notice that before we had m = 1 which means that the particles in 
our system have rest mass 1, but for causal geodesic completeness we need 
to consider any m > 0. Such a geodesic exists on a maximally extended 
interval [0,r + [, and future geodesic completeness means that t + = 00 for all 
such geodesies. 
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The following relations between the variables r, w, L, and p 1 hold: 



Since dt/dr = p° > we can re-parameterize the geodesic by coordinate time 
t G [0,i+[. We remark that since there is no matter in the region r < r' 
the arguments in [2j [23], which apply to any causal geodesic, imply that 

= oo. This will nevertheless be shown below. The arguments in [2J [21] 
are however not sufficient to conclude that r + = co since E and p° may grow 
in time. We now show that E and p°, in the present situation, are bounded 
for any causal geodesic which then implies that r + = co. 

We call the domains [0, r' Q ], [r' Q , Rq], and [Rq, co[ the inner vacuum region, 
the steady state region, and the outer matter region respectively. 

Our strategy is to first show that a geodesic with sufficiently large E, 
for which R < Rq, will at a later time travel outwards for all times. We will 
then apply results from [2j. 

Remark 3: The analysis given below rests in part on the results in the 
previous part of the proof that rapidly outgoing characteristics will continue 
to move outwards. However, since we have good control of the metric and 
the matter in the domain R(t) < Rq + K*t there are alternative approaches 
to the one given below to obtain a bound on E. 

Consider a causal geodesic (R(t), W(t), L). We now re-define the quan- 
tities W*, W- and to adapt these to the causal geodesic we now consider. 
Let 



E 



w 



r 




where we now re-define 



E = E(r, w, L) = \j m 2 + w 2 + L/r 2 . 



W£ := max{l + 




and let W 9 _ > satisfy 



W 9 _e 




2M\ Z1V1 r> In 



Define 
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We have that along a geodesic 



dE 
~dt 



W \ 
At— + ;U r e^ A 



W. 



dW 
~dT 



X t W + fi r e^~ x E - e»- x 



L 

^E 



and 



(3.49) 
(3.50) 

(3.51) 



dt E 

To control solutions of equation (|3.5ip we need to estimate e At_A . We have 
in view of (|3.26j) . 



"° m(t,r) dr 
o r2(l-^g£i; 



R ° M in dr 



Mdr 



, r 1 / 2Mi 
>_41 og _£ + log 1 



r 2 



r 2 (l - ^) y Ro r»(l - ^) 
1 / 2M\ 



Moreover, since the steady state is a given regular solution we have that p 
is uniformly bounded on r 6 [?"o> r o]- This gives 



/2(t 5 r) 



''0 



Ro 



IX, 



> -c 



2A 

4-rrr e pdr 



4-7rr e 2A p (ir 



M 



4irr e 2X p dr 



Ro 



(3.52) 



in analogy with estimate (|3.46p and where ()3.26j) was used to bound e 2A on 
[r ,ro]. Since e~ x = yl — 2m(t, r)/r is strictly positive we have that for 
some constant Cq > 



Define 



e 0i-A)(t,r) > Co for r G [o 
Co 



(3.53) 



P :-- 



m 2 + 1 + L/r'o 



It follows that a geodesic (i?, W, L) with | W| > 1 and R> r' satisfies 

dR 



dt 



> P. 



(3.54) 
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Since the aim is to show that E is bounded we can assume that for some 
to E [0,t+[ 

E(t )>Y ie Y2 + m + ^, (3.55) 
r o 

where Y\ > and where Yi > will be specified below. Note that this 
condition implies that 

\W(to)\ >Y x e Y \ 

We will consider the cases R(to) > Rq; r' < R(to) < Ro; and R(to) < r' , 
and the subcases that W(to) is positive or negative. Consider the first case 
with W(to) positive. In this case we can directly refer to the arguments 
in the previous part of the proof, using the adapted quantities introduced 
above, noticing that G(to) > E{to) > Wl, to conclude that the geodesic 
will travel outwards for all times with R(t) > Ro + n*t. If instead W(to) 
is negative we can assume without loss of generality that the geodesic will 
have R(t) = Ro at some time t. This follows in view of (|3.54p since the 
geodesic must otherwise have |W(t)| < 1 at some time t but then E(t) is 
bounded by y m 2 + 1 + L/ 'Rq and the argument can be restarted at some 
later time to for which (|3.55p holds. Hence we are in the situation of the 
second case which we will treat below. 

Before continuing with the remaining cases we need some auxiliary re- 
sults. The right hand side of (|3.50p depends locally on r except for the 
function /i. Now ji < and since the steady state situated in r' < r < ro 
is a given regular static solution, it follows that there is a constant T > 0, 
such that for \W\ > 1 we have 



dW 



dt 



<T\W\, for r £ [r' , R ]. (3.56) 



Let 

„. T(Ro-r'o, 



P 

We choose Y2 such that Y2 > 2q + 1. We return to the remaining cases. We 
note that in view of (|3.49p . E is constant in the inner vacuum region. Hence, 
in the case where R(to) < r' , we can without loss of generality assume that 
W(to) is positive and R(to) = r' , since if W < the geodesic will continue 
inwards for some time, with constant E, reach a turning point where W = 
if L > 0, and then it will travel outwards, with positive W, and it cannot 
change sign until it hits the steady state region. We remark that if L = 
the geodesic will pass through the origin and the "turning point" will be 
r = where the sign of W changes. Assume now that R(t) < Rq on the 
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time interval I := [to, to + At] where At = (Rq — r' )/P. We then have in 
view of (I330|) that 

log W(t) > log Yl+Y 2 - TAt > log Yi + q + 1, 

for t £ I. This implies in particular, since W(t) > 1, that 

\W(t)\e»- x >p 

on /. In view of (|3.5ip the geodesic has thus necessarily crossed the region 
r' Q < r < Rq within the time interval At. Hence R(t) = Ro with W(t) > 
at some t E I. We can then repeat the arguments in the first case to conclude 
that the geodesic continues to travel outwards. We turn to the case where 
R(to) £ [r' ,Ro]. In the case that W(to) is positive the argument from the 
previous case applies. If W(to) is negative we have by the same argument 
that on the time interval I := [to, to + At] the maximum change for log \ W\ 
is q so that R(t) = r' at some t with log |W(t)| > log Y"i + q + 1. Thus we 
are back in the case where the geodesic is in the interior vacuum region and 
the geodesic will turn and have R(t) = r' Q at some later time t with reversed 
sign on W(t) so that log W(t) >Y\+q + l and we can repeat the previous 
argument. 

To conclude we have reduced the situation so that we only need to con- 
sider geodesies that satisfy R(t) > Ro + K%(t — 1 2 ) for t>t 2 for some t 2 > 0. 
In the case of a timelike geodesic associated to a particle upper bounds on 
G = E + W and H := E — W, and thus on E, are obtained in [2] under the 
assumption that there is no matter in the domain {(t,r) : r < e}. These 
bounds are however time dependent. The arguments in |2] do not depend 
on the rest mass m and are thus unchanged for a causal geodesic. If we 
apply these arguments in the present situation, where the geodesic satisfies 
R(t) > Ro + K%(t - t 2 ), the bounds (4.18) and (4.19) in [2j can in fact be 
seen to be time independent. Indeed, we need to reconsider the estimates 
in [2] leading to these bounds. For simplicity we put t 2 = 0. Let us here 
only consider the bound for G since the bound for H follows analogously. 
We thus reconsider the estimates of (4.14), (4.15) and (4.17) in [2]. For the 
inequality (4.14) in [2] we now get 

f 11 f°° 

\ L R Q + \K<tt\) 

(3.57) 
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Since 

the bound (4.15) in |2] is replaced by 



~We - {r + \kW (3 ' 58) 

Finally, for the term (4.17) in |2] we get 

/ / 8n(p - p)e» +x dtdr < [ /"* — ^1 h2 drds < C(M, R ). 
J Jn Jo JRo+\^\ \ R o + \ K * s \) 

(3.59) 

These estimates turn (4.18) in [2] into the time independent bound 

G(t) < C(M,R ). 

By an analogous argument it follows that also H(t) < C(M, Rq), and there- 
fore E(t) < C(M, Rq). This shows that E is bounded for any causal geodesic. 
Now, since \i is bounded from below it follows that p° and p a are bounded 
and thus t+ = oo, and since dt/dr = p° we have that t + = oo, and the proof 
of the theorem is complete. 

□ 

For the proof of Corollary [1] a slightly different set of initial data is needed. 
The reason is that there is no result in the literature which says that there 
are steady states of the Einstein- Vlasov system for which m/r is arbitrarily 
small everywhere. By numerical simulations, using the code developed in the 
work |llj . we find nevertheless evidence that this is true. It is on the other 
hand known that there are steady states for which 2m /r can be arbitrarily 
close to 8/9, cf. [4J and [3]. 

The set up below is similar to our previous set up with the difference 
that the inner matter which is supported in [j"o,ro] is not given by a steady 
state. 

Let < Tq < ro < r% be given and put M = r*i/2. Let L + > and let 
M out < M be such that 

M - M - < i (3.60) 



'0 2 

Put M in := M — M out . Let R\ > n be such that 



Ri-n< n r ° , (3.61) 
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and define 

Ro := ^(n + Rt). 

Let W* > satisfy 

|W*|>l + ^±. (3.62) 

o o o 

Let / = /in + /out be such that 

supp/ in C [4r ]x]Wi,oo] x [0,1,+], (3.63) 

supp/ out C [R ,Ri}x}W 2 ,^} x [0,L+], (3.64) 
where Wi > and W 2 > satisfy 

_[ 5M in I 5M ^ 

|Wi|e ^ (1 ^f' 2 (l -— -)V2> 3|W2|, (3.65) 



f 3M in , 



'0 

-5M 



) ' 2 «o(i-W); , 2M, 



\W 2 \e ™ (1-— ) 3/2 >3|W,|, (3.66) 

itO 



and such that 



Define 



where 



and 



pro ftii 

M- m = / 4vrr 2 fl n dr, M out = / 4vrr 2 /i ut dr. (3.67) 

=A 3 / 2 B, 



1 + W% + L + /r' 2 



2M in 2M. 

^ := mm 1 — , 1 - — }, 

r Ro 



B :=e ^ 



o 

Theorem 3 Assume that r , r , r\, Rq, R\, L+, M out , W*, Wi, W2 and / are 
given as above, and consider a solution f of the system i2.^) - [27h]) . launched 

o 

by / on its maximal existence interval [0,T[. Then T = 00, and 

supp/(t) C [r / o + |t<|,oo[x[W»,oo[x[0,i+] J (3.68) 
and the resulting spacetime is future causally geodesically complete. 
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Proof: The proof follows to a large extent the previous proof and the set 
up is identical. For a characteristic originating from [Rq, R\], corresponding 
to the outer matter, we consider the same quantity 

G(t)e^\l-^L) (3.69) 

as above. We follow the steps identically until we reach inequality ([3.411) 
where we replace A by A' := {(t, r) : < t < t\, r' < r < oo}. The reason for 
this modification is that there is no matter at r = r' and this fact guarantees 
that the boundary terms which result from the identity (|3.43j) are zero as 
before. The estimate (I3.44p is now replaced by 

oo 



J J re^ +X 2je 2X — drdt < J r 2 e 2X p{0,r) dr 



o 

fRo roo 
< / r 2 e 2X p(0, r) dr + / r 2 e 2X p(0,r) dr 
Jr' JRt) 

1 M- m 1 M 

r' Ho 

Using this estimte in f)3.45j) we get 

4tt J J e » +x pe 2X jdrdt < 4tt J J e» +x 3je 2X jdrdt (3.71) 

< 4vr— / / re^ +X 3je 2X -drdt 
~ Ro J J a r 
3 M- m 3 M 

in place of (|3,45p . The remaining estimates are unchanged and the condition 
(|3.66p then ensures that the characteristic we consider satisfies |iy(ti)| > 
\W m \. 

For a characteristic originating from the interval [rQ,ro] we define the 
time £r such that if the characteristic reaches r = Rq in the time interval 
[0,ii] then this happens at t = tn . Note that on the time interval [0,ti] all 
particles move outwards so the characteristic can only cross r = Rq once. 
On [0, £r ] we consider the quantity 

G(t)e^))(l-^) (3.73) 
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instead of (|3.69[) . The influence from the outer matter can only enter through 
the metric coefficient \i and the actual value of /i plays no role in the proof. 
We can thus follow the steps in the proof above replacing t\ with tn , M 
with M[ Q and Rq with r' until we reach inequality (|3.4ip where we replace 
A by A' := {(i, r) : < t < tn ,r' G < r < oo}. We then again use the 
estimate (|3.72p to replace ()3.45p . The estimate (|3,46p of the boundary term 
is slightly changed in this case since the domain of integration is now [r' , oo]. 
The estimate f|3.46j) now becomes 

e (-M+A)(t,r) A ^ r ) dr+ LOt-m^j^far) -jx t (t, r))dt 

rRo poo 

< / 4irre 2X \j(t Ro ,r)\dr+ / Airr e 2X \j(t Ro ,r)\ dr 

Jr' JR 

< M - + M. (3 74) 

Similarly the estimate (|3.48|) also gives two terms by the same splitting, 
cf. (ET38]) . Hence (ET4T1) becomes 



xG(0)^M°»(l 2Vi ' 



11(0) " 

and inequality f|3.47|) is then replaced by 



G(i iJo )>e 

xG(0)(l-^(l-^)^. 

Here the estimate (|3.48p was modified by splitting the integration domain 
in the subintervals [r-Q,l?o] and [l?o,oo] as above which yields 

/ . , 1 / 2M in \ 1 / 2M\ 
M0,r;)>-log( 1 -^) + -lo g (l-— ). (3.75) 

We use condition f|3.65j) and obtain 

W(t Ro ) > W 2 , 
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by using the arguments in the previous case. 

On the remaining time interval ]t_R ,ii] we again consider the quantity 
G(t)e^ t,R ^ (1 — 2M/R(t)) as for a characteristic originating from [Rq,Ri] 
and repeat the arguments in that situation. Note here that R\ is in this 
case replaced by R\ + t^ but this has no influence on the argument. Thus 
the first statement of Theorem [3] holds for this class of initial data. We next 
consider statement (I3.68H . If we let 

. 2M in 2M 
A : = mm 1 — , 1 - — }, 

r Ro 

we have that e _A( *' r ) > A 1 / 2 . We also have 

eA(tjr) ^ (1 _2Min 1/2 2M 1/2 

in view of inequality ()3.75p . Similarly we get 



analogously to (|3.T4[) . Thus if we define 

^TTW! + L + /r' 2 

it follows that (|3.68p holds. The proof of causal geodesic completeness fol- 
lows the steps in the previous proof but is easier since in this case there is 
no steady state present. 

□ 



4 Proof of Theorem 1 

In this section we prove Theorem [H Corollary [1] and Corollary [2j The 
results of Theorem [2] and Theorem [3] are time reversible in the following 
sense: by taking initial data as specified in Theorem [2] or in Theorem [3] but 
with reversed momenta, disregarding the steady state, so that the outgoing 
particles are ingoing, then global existence to the past holds and (|3.32[) and 
(|3.68p become 

supp/0) C [ai + \a 2 t\, cc[x] - oo,W*} x [0,L + ], (4.1) 
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where a\ and a% are equal to Rq and respectively r' and k*, and spacetime 
is past causally geodesically complete. We denote such initial data by X s 
and X r respectively. We show that there is a sub class of X s and of X r which 
satisfy the conditions in [6J. These conditions guarantee that black holes 
form to the future, cf. also [7] where an additional argument is given to 
match the definition of a black hole in [13]. We recall the set up and the 
conditions on the initial data in [5]. There are two slightly different initial 
data sets in [6] which both guarantee the formation of black holes and for 
our purpose any of these will do. Let us here give the details of Case (i) on 
p.688 in 0. 

Let < r < J"i be given, put M = ri/2, and fix < M out < M such 
that 

2(M - M out ) * 



< n- ( 4 - 2 ) 



Let R\ > n be such that 



and define 



r 9 



Ri-n< ^-A (4.3) 



i? := ^(ri+fli). 

Denote by p the energy density induced by the initial distribution function 
/. It is required that all the matter in the outer region [7*0, oof is initially 
located in the strip [Rq,R\], with M out being the corresponding fraction of 
the ADM mass M, i.e., 

/•oo rRi 

/ 47rr 2 p(r)dr = I 4irr 2 p{r)dr = M out . (4.4) 
Jro J Ro 

Furthermore, the remaining fraction Mj n = M — M out should be initially 
located within the ball of area radius ro, i.e., 

pro 

/ 4irr 2 p(r)dr = M in . (4.5) 
Jo 

If the inner matter is chosen to be a steady state then the solution exists 
for all r G [0,oo[ and for all times but generally this is not required and we 
define the set 

D :={(*, r)€ [0,oo[ 2 | r >7 + (*)}, (4.6) 

where 7 + is an outgoing radial null geodesic originating from r = ro > 0, 
i.e., 

= e 0*-A)(.,^W) +( ) = ro . ( 4. 7) 

as 
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In this set up we use the following expression for the Hawking mass 



m(t,r) = M — / 47T7/ 2 /j(i, rj) drj, 

J r 

instead of (|2.18|) . The analysis in [B] is restricted to D since the dynamics 
of the inner matter is not essential to conclude that a black hole forms. We 
define 

and we require that the parameter W_ < satisfies 

rxn,^) 2 !^! 2 ^, (4.9) 

where 

a ■ J 1 r ° r i- r o\ /, in x 

a = mm<-, , >. (4.10) 

\2' 12fli' 300^1 J v ; 

o 

We impose the General support condition: For all (r, tu, L) G supp / the 
following holds: 

r€]0,r ]U[i2o,i2i], 
and if r £ [i?o ; -Ri] then u; < and also 

3L 

< L < — -m(r/) + rirfiri), rj G [ro, (4-11) 
77 

One of the two classes of initial data specified in [6], which guarantee the 
formation of black holes, can now be given. Let 

{o o 
/| / is regular, satisfies (|4.4p . (|4.5p . the general support condition, 

and for (r,w,L) G supp /with r G [i2o,ili] 3 VZ/r < lj. (4.12) 

Corollary 2.3 in [6] shows that if the inner matter is a steady state then the 
solution exist for all Schwarzschild time. We denote the subset of Ib for 
which the inner matter is a steady state by X S B . The initial data sets I s and 
I r satisfy to a large extent the conditions given above and we are now in a 
position to define the subsets J s C X s and J r C T r which satisfy the claims 
of Theorem [TJ Let 

L+ = min{ rg,r M in }. (4.13) 
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The set J s is defined by 



J s ■= |/ e X s | L+ satisfies (ITOjl and W- satisfies (iOjl j. 

Note that there is no conflict between the conditions ()3.29|) and f|4.9|) by 
taking \W\- sufficiently large. It is clear that J s C X S B . Next we consider 
the initial data set where both the inner and outer matter move rapidly 
inwards. We define 

J r ■= |/= f in + / out e l r | L + satisfies p3| and W 2 satisfies (|PJ) j. 
Since this set is a subset of Z# the proof of the theorem is complete. 

□ 

Proof of Corollary 1: Let e > be given and choose initial data in J T . By 
evolving the data to the past we have by (|3.68p that the matter is supported 
in the domain {(t,r) : r > r' + Thus, by evolving to time t = —T, 

where T is sufficiently large so that 

m(t,r) M 
r r r' + \K*\T 

the claim follows by taking as initial data the solution at t = —T. 

□ 

Proof of Corollary 2: Consider initial data in J s . The statement then 
follows by combining Corollary 2.3 in [6] and Theorem [2j 

□ 
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